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We study the single particle entanglement spectrum in 2D topological insulators which possess 
n-fold rotation symmetry. By defining a series of special choices of subsystems on which the en- 
tanglement is calculated, or real space cuts, we find that the number of protected in-gap states 
for each type of these real space cuts is a quantum number indexing (if any) non-trivial topology 
in these insulators. We explicitly show the number of protected in-gap states is determined by a 
Z"-index, (zi, z n ), where z m is the number of occupied states that transform according to m-th 
one-dimensional representation of the C„ point group. We find that the entanglement spectrum 
contains in-gap states pinned in an interval of entanglement eigenvalues [1/n, 1 — 1/n]. We deter- 
mine the number of such in-gap states for an exhaustive variety of cuts, in terms of the Z m quantum 
numbers. Furthermore, we show that in a homogeneous system, the Z n index can be determined 
through an evaluation of the eigenvalues of point group symmetry operators at all high-symmetry 
points in the Brillouin zone. When disordered n-fold rotationally symmetric systems are considered, 
we find that the number of protected in-gap states is identical to that in the clean limit as long as 
the disorder preserves the underlying point group symmetry and does not close the bulk insulating 
gap. 



The study of novel topological phases of matter has 
become one of the most active fields in condensed matter 
physics^— . Chronologically speaking, the first of these 
phases to be experimentally realized is the integer quan- 
tum Hall (IQH) state, which is a striking departure from 
the traditional theory of conductivity due to its quan- 
tized Hall conductance and chiral edge modest. Shortly 
after its experimental discovery, Thouless demonstrated 
that the special properties of the IQH state came down 
to its non-trivial bandstructure topology^. For the IQH 
state, one can prove that the quantized Hall conduc- 
tance is equivalent to a momentum space integral of the 
Berry curvature, which is a non-zero integer (known as 
the Chern number in topology), multiplied by the con- 
ductance quantum e 2 /h. Beyond this, Haldane demon- 
strated that the IQH state may be further generalized to 
a system which does not require an external overall mag- 
netic flux, yet still possesses a non-zero Chern number 
and quantized Hall conductance. This system is gener- 
ally referred to Chern insulator—. 

Further as systems which preserve time-reversal invari- 
ance (TRI) pose an interesting problem as their Chern 
number vanishes in the presence of non-trivial topology. 
This necessitates the definition of a new quantum number 
capable of distinguishing between trivial and non-trivial 
topological states in systems which possess TRI. In 2Dir— 
and one can define a ./^-number, which here we 

call 7o, as a quantity which is uniquely determined by 
the matrix representations of the time-reversal symme- 
try operator at all time-reversal invariant points within 
the Brillouin zone (BZ). In gapped systems, 70 = 1 
corresponds to topologically non-trivial insulators while 
70 = corresponds to trivial insulators. In real mate- 



rials, TRI topologically non-trivial insulators, or simply 
topological insulators usually have strong spin-orbit in- 
teractions which lead to an inversion of the bulk band 
gap at an odd number of time-reversal invariant points 
within the BZi£^. 

The important difference between Chern insulators 
(characterized by Z-indices) and TRI topological insula- 
tors (characterized by ^-indices) is the absence/presence 
of time-reversal symmetry. In general, we expect that 
the dimensionality and the symmetry of the system will 
determine the available quantum numbers that may be 
used to distinguish the topologically non-trivial states 
from the trivial ones. Motivated by this notion, Schny- 
der et al. classified translationally invariant topological 
insulators in presence/absence of three global symme- 
tries: time-reversal, particlc-hole, and chiral in 2D and 
3 D 26 ' 27 . With translational symmetry, a non-interacting 
insulator maps to a manifold defined on the entire BZ. 
At each k in BZ, the Hamiltonian maps to a projec- 
tor onto the subspace spanned by the occupied Bloch 
states, or an element of the group U(N or b) /U(N occ ) x 
U(N or b — N occ ). The mathematical mappings of var- 
ious insulators have been previously studied and it is 
found that they may be classified by their specific ho- 
motopy groups^. To be more concrete, in 2D the ho- 
motopy group for systems with only translational sym- 
metry is n 2 (U(N orb )/U(N occ ) x U(N orb - A^ occ )) = Z, 
which indicates that topological phases are indexed by 
integers. This picture is consistent with the previously 
understood conclusion that 2D systems may be charac- 
terized by a Chern number. In the same manner, TRI 2D 
and 3D insulators are also classified by the correspond- 
ing homotopy group, Z2, of continuous mappings from 
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the d = 2, 3-dimensional BZ to a target space of pro- 
jectors to occupied states under the constraint of TRI. 
Besides the bulk topological indices defined in insulators 
with translational invariance and with disorder, a bulk- 
edge correspondence exists which guarantees presence of 
gapless modes on the boundary between two insulators 
having different values of Z or Z% inde x 25 ! 29 ' 30 . 

The three symmetries considered above are local sym- 
metries that do not involve spatial degrees of freedom. 
In condensed matter systems, we are faced with many 
types of symmetry operations in real space, besides 
time-reversal, charge conjugation and chiral symmetries. 
Amongst these symmetries, crystallographic point group 
symmetries (PGS) are well-known to universally exist 
within solids and may thus provide a natural avenue 
through which to continue the expansion of the family of 
materials which harbor topologically non-trivial states. 
Along this line, Fu 3 ^ first studied spinless 3D insulators 
exhibiting gapless surface modes with quadratic band dis- 
persion, the existence of which is protected by C/±$ rota- 
tion PGS and TRI. Furthermore, within the context of 
cold atoms, Sun et al. discussed a quadratic band cross- 
ing point in 2D that is protected by C4 PGS and TRI. 
These examples provide a tantalizing glimpse into the 
possibilities which exist for finding new non-trivial states 
of matter when PGS is combined with TRI. An equally 
interesting question one may ask is: can PGS bring about 
new topological phases, without the presence of TRI and 
how shall we identify the m 17 ' 19 ? Without the presence 
of TRI, we are no longer a priori guaranteed to have any 
of the topological characteristics which are hallmarks of 
this symmetry such as: quantized magneto-electric coef- 
ficient or boundary /surface modes thereby necessitating 
the definition of a new topological quantum number. 

To this end, the study of insulators with inversion 
symmetry- 3 - 2 - - — , which is the simplest yet non-trivial 
PGS, shows that in these insulators a new topological 
number may indeed be defined which is independent of 
the Z- or Zi numbers previously introduced. This num- 
ber can be used to characterize an inversion invariant 
insulator which contains a null result for both the Chern 
number (Z-number) and ^-number, but is still topo- 
logically distinct from a trivial insulator. The new quan- 
tum number labeling the new topological insulators is the 
number of protected in-gap states in the single particle 
entanglement spectrum. 

The entanglement spectrum has been successfully ap- 
plied to identify topological orders in a variety of con- 
densed matter systems such as fractional quantum Hall 
systems 3 - 5 -, spin chains 3 ^, Chern insulators 3 - 7 -, and topo- 
logical band insulator s 30 ' 38 " — . Unlike the energy spec- 
trum, the quantum entanglement spectrum solely de- 
pends on the many-body ground state, and is, for this 
reason, more suitable for identifying any non-trivial 
topology in the ground state. In noninteracting systems, 
it can be shown that to obtain the many-body entan- 
glement spectrum, only the single particle entanglement 
spectrum is needed, defined as all eigenvalues of the one- 



body reduced density matri x ' . 

a, (A) = { c \c 3 ), (1) 

where i,j £ A and A is a sub-system of the whole system 
L. In principle, all 'energies' in the entanglement spec- 
trum are between zero and unity corresponding to states 
which arc cither localized outside of or inside the sub- 
system, respectively. In the case of inversion symmetric 
topological insulators 3 - 2 -, if A is chosen as the left/right 
half of the L, then there can be several eigenstates of 
C(A) with degenerate eigenvalue at exact 1/2, in an in- 
version invariant insulator. These degenerate eigenvalues 
correspond to states that are topologically protected in 
the sense that they always stay at pinned 1 /2 so long as 
the inversion symmetry is preserved and the bulk gap re- 
mains open. These states are referred to as the protected 
1/2-in-gap states and their number defines a new topo- 
logical number to be used to differentiate between trivial 
and non-trivial inversion invariant insulators. 

In this work, we focus on insulators with another type 
of PGS, n-fold rotation symmetry where n = 2,3,4,6 
by lattice restriction. We show that for proper choices 
of sub-system A in PGS invariant insulators, there ex- 
ist protected in-gap states the number of which defines a 
new topological invariant. While these states are analo- 
gous to the 1/2-in-gap states previously seen in inversion 
invariant insulators, two key differences are noted: first, 
instead of taking A as the left /right half of L, there exists 
more than one type of symmetric cuts for a C n invariant 
insulator and, in general, each has its own corresponding 
number of in-gap states. Second, the protected in-gap 
states now are not necessarily located at exact 1/2, but 
within a region around 1/2, prevented by symmetry from 
arbitrarily approaching two ends of the single particle 
entanglement spectrum (zero and unity). We calculate 
the number of protected in-gap states for each type of 
symmetric cuts in a C n invariant insulator in terms of 
m-fold rotation eigenvalues at high-symmetry points in 
the BZ, where to is a factor of n. Further, we show 
that all the numbers of protected in-gap states can be 
expressed in terms of a i? n -index (z\, Z2, z n ), giving a 
^"-classification of these insulators. The physical rea- 
son underlying the Z n classification is very simple: if the 
single particle Hamiltonian is 71-fold invariant, then ev- 
ery single particle wavefunction belongs to a certain ID 
representation of C n group. The number of electrons in 
each representation below the Fermi level then naturally 
gives the index (zi, Z2, z n ). Finally, we will show that 
if the translational symmetry is broken by introducing 
disorder potential that is also C„ invariant, the i?™-index 
remains unchanged and as do the associated numbers of 
protected in-gap states. 

The paper is organized as follows. In Sec[H we intro- 
duce concepts related to our examination of PGS insula- 
tors in 2D referred to throughout the paper: In SecOjA), 
we introduce crystallographic point groups; in in Scc|T}3, 
we discuss a general definition of non-interacting topo- 
logically non-trivial insulator, and in Scc|T|C), a brief in- 
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traduction to single particle entanglement spectrum and 
its relation to entanglement entropy. In Sec|H]we revisit 
the entanglement spectrum in a ID inversion invariant 
system. In this revisit, we analytically prove the rela- 
tion between the number of protected in-gap states and 
the inversion eigenvalues at k = and k = ir. The em- 
phasis here is the relation between the number of 1/2-in- 
gap states and a Z 2 -index (zi,Z2), where z\ (22) is the 
number of occupied Bloch states with odd (even) parity. 
Sec lIIII we perform a simple extension of the previous ID 
results to 2D and 3D inversion invariant insulators, where 
we relate the number of 1/2-in-gap states to the parity of 
Chcrn number and, when considering systems with TRI, 
to the Zi number. In Sec lIVl we present our main result. 
We outline the analytic proof by which we relate the Z n 
index to the number of protected in-gap states in a C n 
invariant insulator with every type of possible symmetric 
cut, with some specific details relegated to the various 
Appendices. In Sec[V] we briefly discuss 3D PGS invari- 
ant insulators in SecEjA); we discuss in Sec|VjB) the 
effect of weak interaction, which reduces the Z"-indcx 
down to a Z„-index as single particle states lose their 
meaning while the many-body wavefunction of any non- 
degenerate insulating state is still a ID representation 
and can be associated with some Z n number. 



I. PRELIMINARIES 

Before proceeding with the relevant background ma- 
terial, we briefly clarify the notations we use throughout 
the paper. When we use the notation C n , we are referring 
to either the point group C„, or to the n-fold rotation, 
the choice of which should be clear given the context. C n 
refers an operator in the Hilbert space corresponding to 
the n-fold rotation. C n denotes the matrix representa- 
tion of the operator C n in a basis, either the orbital basis 
or the occupied-band basis, which should be discernible 
given the context of its usage. Similar definitions and 
distinctions apply to other operators, e.g., M Zl M z , M z 
refer to mirror symmetry z — > —z, its operator in Hilbert 
space, and a matrix representation of the operator, re- 
spectively. Finally, when used as indices, the Greek let- 
ters will always span the orbitals, while Roman letters 
denote either the bands or spatial directions (x, y, z or 
k x ,k y ,k z ). In all equations, repeated indices are, unless 
otherwise noted, automatically summed for the follow- 
ing types of summations: summation over all sites in the 
system, summation over all sites in the subsystem on 
which the entanglement spectrum is calculated, summa- 
tion over all orbitals, summation over all bands (occupied 
and unoccupied). 

A. 2D Crystallographic Point Groups 

A point group is a group of all symmetry operations 
(an atom, a molecule or a lattice) that leave at least one 



point fixed in space. In nature, there exist an infinite 
number of point groups. On the other hand, crystallo- 
graphic point groups are point groups that are consis- 
tent with lattice translational symmetry. The crystal- 
lographic restriction theorem states that lattice trans- 
lational symmetry dictates that n = 1,2,3,4,6 for any 
n-fold rotation axis and effectively limits the total num- 
ber of crystallographic point groups. In 2D there are 
10 different point groups allowed which are comprised 
of five cyclic groups and five dihedral groups. A cyclic 
point group C n consists of all powers of an 71-fold rotation 
while a dihedral point group D„ contains C„ as a sub- 
group and an additional two-fold in-plane rotation axis. 
In this paper we will focus on C„ PGS. 

With the point groups defined, we now consider the 
representations of point groups in a single particle Hilbert 
space. In the representation, every symmetry operation 
corresponds to an operator in the Hilbert space. As we 
will discuss both spinless and spinful fermions, it is nec- 
essary at this point to introduce double point groups, 
G D . Mathematically speaking, G D = {E,E}®G, where 
E and E are the identify operator and a rotation of an- 
gle 2ir about any axis in space, which satisfies E 2 = E, 
and Eg = gE for g E G. This enlarged group has 
physical meanings as one considers its representations, 
which are classified into two types: in the first type 
D(E) = D(E) = I, where I is the identity matrix; 
in the second type D(E) = —D(E) = 1. Represen- 
tations of the first type are referred to as single- valued 
and representations of the second type are called double- 
valued. A degenerate subspace, in which every state 
has the same energy, of a single fcrmion must form a 
single- valued/double- valued representation of the under- 
lying PGS of the Hamiltonian if the particles are spin- 
less/spinful, as we know that a spinful fermion takes an 
extra minus sign as it is rotated by 2n. For C„ PGS, as 
one applies the general statement to Hilbert space rep- 
resentation of C n PGS, one can obtain (7™ = ( — 1) F I, 
and if in the orbital/band space representations, simi- 
larly, CI = {-1) F I, where F = (F = 1) for spinless 
(spinful) fermions. 

Given the Hilbert space representation R of a point 
group operator R, we now look at the sufficient and neces- 
sary conditions for a single particle tight-binding Hamil- 
tonian H to be invariant under this operator. We assume 
that in one unit cell at R = niai + n2^2 + n^Si^ there 
arc s orbitals at r = R + d, where d is the offset vector 
of the atom from the nearest lattice point. Note that we 
have ignored all unit cell structures, a simplification that 
spares us from discussing non-symmorphic space groups, 
which contain operations combined from member of a 
point group and translations by fractions of lattice vec- 
tors. Now the PGS operator R sends an electron at r to 
Rr = RH + Rd. If R is a symmetry, there must be an- 
other atom at Rr, i.e., R' = RR + Rd — d is also a lattice 
point. Writing the above relation in second quantization, 
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wc have 



to p z . In this case we have in the basis of (p x ,Py,Pz) 



Rc a (R)R- 



(2) 



0-10 
724 = | 1 
1 



(3) 



In this equation, lZ a p describes a possible rotation of an 
orbital under R. For example, suppose we are considering 
Px,j/,2-orbitals under R = C4 a four-fold rotation about 

z-axis, then obviously, p x is sent to p y , p y to — p x and p z from Eq.@, we obtain: 



If we adopt the following definitions of Fourier transform 
of c Q (R + d) 

c Q (k) = ^^ CQ (R)cxp(-zk-R), (4) 



RcaO^R' 1 = -/=Y1 K a pcp (R') exp(-ik • R) (5) 
= -== K^cpCR') exp[-ik • (R- 1 ^' + iT M - d)] 
= -L V 72 Q( 9C /3 (R') exp[-i(i?k) • (R + d - i?d)] 

^ N R 

= 7e Q/3C ^(i?k) e lk ( d -« d ) 

= 72 Q/ 3(k) C)3 (i?k). 



It is easy to verify that 72(k) is a unitary matrix us- 
ing that R is a unitary symmetry operation. A tight- 
binding model Hamiltonian with translational symmetry 
H = X)k^ ct /3(l <; )c3i(k)c/3(k) transforms under R as 

RHR- 1 = (72.(k)^(k)72^ 1 (k)) ct(3 c Q (i?k)c (3 (i?k)(6) 

k 

If R is a symmetry, we have RHR^ 1 = H therefore 

£(k)W(k)£ _1 (k) = «(iZk). (7) 

Using ft(k) = 72e ik '( d - fld ), we have 

KH^K- 1 =H(Rk). (8) 

Eq. (J8J! determines the sufficient and necessary conditions 
for a single-particle tight-binding Hamiltonian to be in- 
variant under a point group symmetry R. For a C n in- 
variant 2D system, the above equation simplifies to 

C n H{k x ,k v )C- x = (9) 

, . 2TT . 27T . 27T /27T.. 

rt(fc x cos( — ) — k v sm( — J, fe x sm( — ) + k y cos( — )), 
n n n n 

where C n is the transformation matrix in the orbital basis 
for the 71-fold rotation about z-axis. 



B. Generalized definition of a topological insulator 

In this paper, we define a topological insulator as a sys- 
tem the many-body ground state of which cannot be adi- 
abatically deformed into the ground state in the atomic 



limit^i^. The atomic limit is the limit of zero coupling 
between orbitals on different atoms but with all symme- 
tries of the system preserved. Every system in this limit 
is considered as a trivial insulator. Clearly then, insu- 
lators that cannot be continuously tuned to the atomic 
limit while preserving the gap and all symmetries are 
nontrivial topological insulators according to this defi- 
nition. This definition enables us to use any quantized 
quantity depending on the wavefunction to differentiate 
topologically trivial and non-trivial insulators. This is 
because a quantized quantity only changes when gap is 
closed, thus if the value of this quantity in an insulator 
is different than that in its atomic limit, it must be a 
topologically non-trivial insulator, or simply, a topologi- 
cal insulator. Such a quantity is called a topological in- 
dex, the choice of which is determined completely by the 
dimensionality and the given symmetries of a Hamilto- 
nian. This definition comprises more types of non-trivial 
topologies than those that may be labeled by a non-zero 
Z or Z2 index. Furthermore, it should be noted that a 
good quantum number that denotes a topological state 
is not necessarily a physical observable such as Hall con- 
ductance. Such a case has been noted previously in non- 
trivial topological insulators possessing inversion symme- 
try. These systems have topologically protected in-gap 
states at exactly 1/2 in the entanglement spectrum and 
thus they cannot be removed by adiabatic deformations 
of the Hamiltonian^. Nevertheless, while the number of 
in-gap states is a topological index, it is not related to 
any yet known experimental observable. 
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C. Single Particle Entanglement Spectrum 

We have defined the general single particle entangle- 
ment spectrum in Eq.([T]), and here we restrict the def- 
inition to the entanglement spectrum associated with a 
special type of subsystem: real space cuts. To begin, the 
correlation matrix for a real space cut is given by 

C MJ/3 = (ct (r(i)) C/3 (r(j))>, (10) 

where r(i) denotes a unit cell within a predefined subset 
of the system, called A, and a, f3 denote orbitals in a unit 
cell. Suppose in the selected subset one has Na unit cells 
and each unit cell has N or b orbitals, the entanglement is 
then given by a Na * N or i, x Na * N or i, matrix, where 
Na is the number of sites in the real space cut A. The 
set of eigenvalues of the matrix is defined as the single- 
particle entanglement spectrum or simply entanglement 
spectrum associated with subset A. The definition is a 
special one as it confines the c-opcrators within a pre- 
defined set of sites in real space. In general, there are 
other types of cuts such as momentum cuts or orbital 
cuts 3 -^, which confine the c-operators within a range of 
momentum or a set of orbitals (e.g., Landau orbitals). 
From definition, it is easy to show that all eigenvalues 
of C, i.e., all 'energies' in the spectrum, are within the 
range [0, 1]. A cut may be chosen as the whole system, L, 
and one can prove that in this case the correlation ma- 
trix C(L) must be a projector matrix, whose eigenvalue 
is cither zero or unity 3 - 2 -. 

The reason we are only interested in real space cuts in 
this paper is that we will focus on band insulators, whose 
constituent single particle wavefunctions are completely 
local in momentum space with k as a good quantum num- 
ber. Therefore all eigenvalues in a momentum cut must 
be either zero or unity, a situation usually referred to as 
trivial. On the other hand, a real space cut has a po- 
tential to distinguish insulators that can be adiabatically 
tuned into the atomic limit and those that cannot be. To 
better understand this point, let us consider an insulator 
in the atomic limit. In this limit all single particle wave- 
functions are localized, and therefore when a real space 
cut is made, a particle must be either totally inside or 
totally outside the subsystem defined by the cut. This 
forces the eigenvalue of the correlation matrix to be ei- 
ther 1 if the particle is located within the cut or if the 
particle is outside of the cut. If an insulator cannot be 
adiabatically tuned to the atomic limit, then within the 
entanglement spectrum there must be eigenvalues that are 
topologically prevented from approaching arbitrary vicini- 
ties ofl or 0. Wc may gain additional insight into the en- 
tanglement spectrum by considering a 2D IQH state, or a 
Chern insulator, that has C (=Chern number) protected 
extended states and, when a real space cut is made, there 
must be at least the same number of states that are nei- 
ther inside nor outside the subsystem defined by the cut. 
These extended topological states present themselves as 
protected 'in-gap', neither nor 1, states in the entangle- 
ment spectrum. This means that existence of protected 



in-gap states in the entanglement spectrum implies exis- 
tence of intrinsically delocalized states. 

The single particle entanglement spectrum is closely 
related to the entanglement entrop y 42 ' 44 " — by 

E q = ~ ln&, + (1 - £„) ln(l - £„)], (11) 

n 

where is the nth eigenvalue in the spectrum. Entan- 
glement entropy is known to quantify non-trivial topol- 
ogy in insulators: If the entanglement entropy cannot be 
adiabatically tuned to zero, the insulator is topologically 
non-trivial. When there are protected in-gap states in the 
entanglement spectrum, then from Eq. (|ll[) the entangle- 
ment entropy always remain finite. Therefore, presence 
of protected in-gap states indicates non-trivial topology 
in an insulator. 



II. ENTANGLEMENT SPECTRUM OF 
INVERSION INVARIANT TOPOLOGICAL 
INSULATORS IN ID: A REVISIT 

A. Decomposition of the correlation matrix into 
projectors 

Let us begin our study of the entanglement spectrum 
of point group symmetric topological insulators by ex- 
amining the inversion invariant syste m 32 ' 34 in ID. The 
system is assumed to have N sites where N is an even 
number. The inversion center is at the midway of the 
iV/2-th site and (N/2 + l)-th site. For now, we assume 
that the system is translationally invariant, hence from 
Eq.© we have [P,H{0)} = [V,H(ir)] = 0, by which we 
can label all eigenstates at and tt with their eigenvalues 
under V (±1). k = and k = tt are called the invariant 
fc-points under inversion, denoted by fcj n „. For a reason 
that is detailed in Appendix we now use, instead of 
eigenvalues of V, eigenvalues of V{k) = 'p e lk i d - pd ) to 
label the eigenstates at k — 0, tt. For inversion invariant 
system, the symmetry fixed point is in the middle of two 
sites so d = 1/2, and d - Pd = 1/2 - (-1/2) = 1. From 
this we have V(k) = e lk V. 

Now we separate the left-half of the system, from the 
1st to the A/2th site, as the subsystem A and calculate 
the entanglement spectrum {£„|A}. There can be eigen- 
values in the spectrum lying exactly at £ ra = 1/2, and the 
number of these eigenvalues is given by 

Ni/2 = I ^2(n + (k inv ) - n-(k lnv ))\, (12) 

k inv 

where n + (ki nv ) and n~(ki nv ) represent the number of 
occupied states with +1 and —1 eigenvalue oiV(ki nv ) of 
the mth band, respectively. If there is an impurity po- 
tential, it will generally split the 1/2-in-gap degeneracy; 
however, if the impurity potential is also inversion sym- 
metric, having equal magnitude at the i-th site and the 
(N + 1 — i)-th site, the mid-gap states are preserved. 
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(c) (d) 

FIG. 1. Examples of inversion symmetric cut in a ID system 
with six sites. All sites in red define the real space cut on 
which the entanglement spectrum is calculated. 

In fact, we desire to make the statement more general 
by showing that the same result holds for all 'inversion 
symmetric cuts', of which the above left-half cut is a spe- 
cial case. The subset which defines an inversion symmet- 
ric cut does not have to be the left or the right half, but 
should satisfy the following conditions: First, the subset 
A has exactly N/2 sites, and second, if site i is in A, site 
N + 1 — i must not be in A. Such a subset A will be re- 
ferred to as an inversion symmetric cut. A set containing 
all sites left to the inversion center is a symmetric cut, 
and several other possible symmetric cuts are plotted in 
Figffl Below, we will use r.i to denote the position of the 
i-th site in the cut for i = 1, N/2. When the sum is 
over i, it means summing over rf, when the sum is over 
r, it means summing over all sites in L. 

While the correlation matrix C in Eq. (fTOj) is not a pro- 
jector, if A is an inversion symmetric cut in an inversion 
invariant insulator, C can be expressed as the linear av- 
erage of two projectors D and D. To show this, we define 
new operators: 

d ia = (c a {r t ) + Pc a {n)P)/V2, (13) 

dia = {c a ( ri )~ Pc a {n)P)/V2, 

where i = l,...,N/2 and a is the orbital index (includ- 
ing spin). These are the symmetric and antisymmetric 
combinations of the annihilation operators. The single 
particle entanglement matrix for subsystem A is: 

c^p = (cUn)c e (n)) (14) 

= ^(4+4)(^+40>- 

Then we point out a simple fact: 

(dtJjff) = (d] a d lfj ) = 0. (15) 

To prove this one only needs to notice that the ground 
state wavefunction must be an eigenstate of P, and that 
PdiaP = d ia while Pd ia P = —d ia . Using Eq. ([T5l) . we 
have 

C ta ,j/3 = ^{Dta,jl3 + £>ia,jp), (16) 

where 

D ia ,jp = (d\ a d jP ), (17) 

Aa,j/3 = (d\ a djp). 



Next we prove that D and D are projector matrices, 
i.e., D 2 = D and D 2 = D. In an inversion invariant 
system, all single particle eigenstates can be divided into 
two subspaces with even and odd parity respectively. We 
use %j) n to denote the annihilation operator of an even 
parity eigenstate and ip n to denote that of an odd parity 
eigenstate, i.e., Pi^ n P = 4>n and P"ip n P = —4> n - Here 
n is not a band label as before, because we do not have 
translational symmetry, but simply an index counting all 
single particle states. In a system with TV sites and N or b 
orbitals per site, there are NN or b/2 states with even and 
odd parity respectively, so n ranges from 1 to NN or b/2 
in general. Every creation/annihilation operator can be 
expressed in this basis and we have the unitary transform 

Then from Eq. (|13[) we have d ra = \2U ra>n tp n . On the 
other hand we have 

ipn = U^ ra c ra , (19) 

tfj n = P^J n P (20) 

— r/t Pc P 

= UnrfjPpaCN+l-ra- 

Equating right hand sides of the last two equations, we 
have 

U ra ,n = 'Pal3UN+l-rf3,n- (21) 

Using these relations we can prove that D is a projector 
by calculating D 2 directly. In fact, 

D j p, ia = (d] /3 d ia ) = 2 Ul >r . U ria , m . (22) 

mGoec. 

And we have 

Djp >ia Dt y> jp = 4 ^2 U r ia ,mU mir .pU r:j p >n Ul triV 

m,n£occ. 

(23) 

where the summation is over the sites belonging to the 
A-subsystem. 

Care should be taken to note that r .pU rj p tn ^ 6 mn \ 
This is because by summing j we only sum half of the 
system. When inserting Eq. (f2Tj) we obtain: 

U li, rj p U r 3 0.n = Ul lN+1 _ r]Ct 'P' al3 'Pp^U N+ l- r ^, n (24) 

= Um^+i-rjpUN+i-TjP^ [Do not sum over j]. 

This leads to 

Substituting Eq.® into Eg. ([23]), we see D 2 = D, or D 
is a projector. 

Using Eas. (|22l25|) . it is easy to see that every row in 
U rjan is an eigenvector of D. If ip n annihilates an oc- 
cupied state, this vector has eigenvalue 1 and if ip n an- 
nihilates an unoccupied state, it has eigenvalue 0. So 
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the dimension of D, or the number of unity eigenvalues 
is: dim(D) = Dim(^f occ ), or the number of occupied 
cigenstatcs with even parity. Similarly, it can be shown 
that dim(D) = Dim(^ occ _), or the number of occupied 
eigenstates with odd parity. 

If dim(D) > dim(D), there must be at least dim(D) — 
dim(D) common eigenstates of D and D with unity and 
zero eigenvalues, respectively. These states are there- 
fore eigenstates of C with eigenvalue exactly at 1/2, 
i.e., 1/2-in-gap states. And if dim(D) < dim(D), there 
must be dim(D) — dim(D) 1/2-in-gap states. In gen- 
eral, we have N\n = \dim(D) — dim(D)\. The two in- 
tegers dim(D) = Dira{^ occ ) and dim(D) = Dim(^> occ ) 
are the numbers of occupied eigenstates with even and 
an equal number of states with odd parities, or in other 
words, they are eigenstates in the two ID representa- 
tions of the inversion point group. These two numbers, 
(zi,z 2 ) = (Dim(ty occ ),Dim(^ occ )), give a Z 2 classifi- 
cation of inversion invariant insulators. The sum of its 
two components z\ + z-i is simply the total number of 
fermions, and their difference \z\ — Z2I is exactly the num- 
ber of protected in-gap states in the entanglement spec- 
trum, i.e., 

N 1J2 = \z x -z 2 \. (26) 

In a trivial insulator, one has z\ = z 2 , while any pair 
{z\,z 2 ) with z\ ^ z 2 denotes a topologically non-trivial 
inversion invariant insulator. Moreover, it should be 
noted that the proof presented here, up to this point, 
does not require translational symmetry. 

B. Translational invariance and number of in-gap 
states 

Now we may proceed to calculate Dim{ty occ ) and 
Dim(^> occ ) by relating them to the inversion eigenval- 
ues at high symmetry fc-points and n in presence of 
translational symmetry. Suppose 7 m (fc) is the annihila- 
tion operator of the Bloch wavefunction at k on the mth 
band, and define symmetry adapted operators 

lM*0 = (7m(*) + Plm(k)P)/V2, (27) 
$m(fc) = (7m(*) - Pj m (k)P)/V2, 

and they satisfy Ptp m (k)P = %j) m {k), P^ m {k)P = 
—ip m (k). Due to inversion symmetry, the single particle 
states \i/}(k)) = ^ f (k)\0) and |^(fc)) = ■0 t (fc)|O), in which 
|0) refers to the vacuum state and not the Fermi sea, are 
the eigenstates of H. The two sets together, for m e occ, 
make a basis of the filled bands in which every basic vec- 
tor has a certain parity. The subspaces ^itP and ^icP 
are given by: * occ _= {\ip m (k))\m <E occ.,fc > 0, ip m (k) ^ 
0} and * occ = {\i/) m (k))\m G occ.,k > 0,t/) m (k) ^ 0}. 
We observe that we have a rather strange looking con- 
straint ip m (k) 7^ and ip m (k) ^ 0. This is because at 
k = 0,7r, either ip m {ki nv ) = or 4> m (ki nv ) = 0, and ei- 
ther the state \i/j m (k)) or \4> m (k)) becomes a null state. 



|V(fc)) + m-k)) G *, 




— 7T —k k 7T 

|V(0/7T)> G * QCC ifS(0/7T) = l, 

|V(0/tt)) e * occ iffi(0/7r) = -l, 

FIG. 2. A schematic of a single band in a ID inversion sym- 
metric insulator. From the figure we clearly see that at any 
k = we have two equal energy states at k and — k, while at 
k — 0, 7T, there is only one state. The two degenerate states 
at ±fc can be recombined into two states with even and odd 
parities respectively, while the singlet state at or it is either 
even or odd. 



To see this, we introduce the sewing matrix 

Bmn{k) = (u m (-k)\P\u n (k)) (m,n G occ). (28) 

At k — and k — n, Hamiltonian commutes with inver- 
sion, and, therefore, the sewing matrix must be diagonal 
with diagonal elements being either +1 or — 1. Apply- 
ing the sewing matrix by utilizing a property proven in 
Appendix [Bj we have 

P 7m (fc)P = B nm7n (-fc). (29) 

Inserting the above equation into Eq. (|27|) and take k = 
0, 7T, we obtain 

^m{kinv) — (1 ~t~ Bmm{.kinv))~{m{,kinv) / 'v2 3 (*^0) 
^m{kinv) (1 $mm{kinv))^/rn(kinv) / '"n/2- 

It is then obvious that if B mm (ki nv ) = 1, then 
i>m{kinv) = 0, otherwise ip m (ki nv ) = 0. This is the 
mathematical description of the rather simple fact that 
for each at k ^= 0,7r, there is its counterpart at — fc, from 
which we can construct two states with even and odd 
parities respectively, while at k = 0,7r we cannot do this 
and either the even or the odd state must be a null state. 
The counting of states in ^ occ and "J occ can be done on 
each band separately and then we add the numbers to- 
gether. At k 7^ 0, 7T, tp m (k) and ipm(k) appear in pairs 
so they contribute equally to Dim(^ occ ) and Dim{^> occ ) 
(see Fig(2]for an example). The number of states in \l/ occ 
and ^occ contributed by this part is (N — 2)/2 for each. 
Counting the states at inversion invariant points needs 
information of the sewing matrix. 

To be more concrete, let us examine a system with 
only one band filled. Suppose we have i?(0) = B(ir) = 1, 
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using Eq. ([30)) . then 

V>(0) = = v/2 7 (0), (31) 
^(0) = 0, 

i>(ir) = 0. 

Both at k = and k = tt, there is one more state in ^ occ 
than in ^ occ - From this we see that z\ = Dim(fy occ ) = 
(N — 2)/2 + 2 and z 2 = Dim(V occ ) = (N - 2)/2. Ac- 
cording to previous discussion, we have \z\ — z<^\ = 2 
1/2-in-gap states. One can repeat the discussion for 
5(0) = B{n) = -1 and hnd z x = Dim(^ occ ) = (N-2)/2 
and z 2 - Dim(^ occ ) = (N — 2)/2 + 2. There are also 
two 1/2-in-gap states for this case. If B(0) = —B(ir), wc 
have z\ = z 2 = N/2, which means no protected 1/2-in- 
gap states. 

For more than one filled band, we repeat the process 
for each filled band, and find that z\ = Dim(^> occ ) is the 
total number of +1 in the diagonal of £>(0) and B(ir), 
and z 2 = Dim(^> occ ) is the total number of —1 in the 
diagonal of B(0) and B(ir): 

z\ = n + (0) + n + (7r) = y] n + (k inv ) (32) 
z 2 = 7Z_ (0) + n_(7r) = ^2 n-(k inv ). 

Substituting these expressions into Eq.(|26j). we recover 
Eg. (1121) . (The sewing matrix used here is different from 
the one used in Ref . [32| due to different choice of inversion 
center.) 

C. Protection of in-gap states against inversion 
symmetric and weak disorder 

We would like to understand how these new quan- 
tum numbers, z\ and z 2 will change when we consider 
a disordered system. Since our previous discussion in 
SecllIjA) did not rely on translational invariance the an- 
swer will not change so far as the disorders are inver- 
sion symmetric. We begin by assuming a form for the 
on-site potential which preserves the inversion symmetry 
V = E ar V{r)c\(r)c a (r) satisfying V(r) =V(N+l-r), 
then we have PVP — V . Therefore, 

(Vw|F|V>„> = 0. (33) 

This entails that if V is tuned adiabatically, an even par- 
ity state \iprn) G ^ will always remain in the even par- 
ity subspace Vf, so the total number of occupied states 
of even parity, Zi, does not change during the process. 
Similar argument gives that z 2 also remains unchanged 
during this adiabatic process. The many-body ground 
state for H + V will be the same as the state evolved 
from the many-body ground state for H, if there is no 
level crossing in the adiabatic process. Therefore, we 



conclude that, as far as \V\ is small compared with the 
bulk gap, the numbers of occupied states z\ and z 2 in 
even and odd subspaces are the same as z\ and z 2 in the 
homogeneous system. Thus, the number of 1/2-in-gap 
states, \z\ — z 2 \, is unchanged and the 1/2-in-gap states 
are robust against any inversion symmetric disorder, and 
do not depend on translational symmetry in general. 

Before ending our revisit to the ID inversion invariant 
insulators, we must emphasize the following two points 
that give us insights for extending our work on inversion 
invariant insulators to C n invariant insulators. First, we 
want to choose a subsystem that is conjugal to its com- 
plementary subsystem with respect to the point group 
symmetry, because in this case, the correlation matrix 
decomposes into a linear average of projectors. Second, 
we can define an integer for each ID representation of 
the point group, which is the total number of occupied 
states that transform according to that representation. 
The sum of these integers gives the total number of elec- 
trons, while their difference is related to the protected 
in-gap states in the entanglement spectrum. 



III. ENTANGLEMENT SPECTRUM OF 
INVERSION INVARIANT TOPOLOGICAL 
INSULATORS IN 2D AND 3D 

The conclusion we arrived at for the ID insulators with 
inversion symmetry can be easily generalized to higher 
dimensions. In these insulators, one can also define the 
Z 2 index as the number of occupied states with even and 
odd parities: (zi,z 2 ) = (Dim(^ occ ) , Dim{^! occ )) . The 
sum of its two components z\ + z 2 gives the total number 
of electrons and their difference, \z\ — z 2 \, is the number of 
1/2-in-gap states in the entanglement spectrum for any 
inversion symmetric cut. Similar to the symmetric cut in 
ID, one defines a symmetric cut in 2D (3D) as a proper 
subset of sites called A that has exactly one half of all 
sites, and if site r £ A, then PrGA. The detailed proof 
of N1/2 = \zi — z 2 1 is omitted because it takes exactly the 
same steps, which we simply summarize here as: 1. Write 
the correlation matrix C as the linear average of two 
projector matrices D and D. 2. Show that the number 
of 1/2-in-gap states is the difference between the number 
of unity eigenvalues of D and D. 3. Show that dim(D) = 
Dim(ty occ ) = z\ and dim(D) = Dim{^> occ ) = z 2 . 

We can easily count Dim{^ occ ) and Dim(ty occ ) in 
presence of translational invariance. For each pair of 
(k, — k) where k ^ ki nl) , there are for each band ex- 
actly one parity odd and one parity even eigenstate, thus 
contributing equally to Dira(^ occ ) and Dim(^f occ ) . The 
contribution from all these pairs is (N — A)N occ /2 in 2D 
and (N — 8)N occ /2 in 3D, as there are four and eight 
k inv 's in 2D and 3D respectively. In the following, we will 
neglect the contribution from all of the aforementioned 
pairs in the Z 2 index, as they are completely determined 
by the total number of sites and the total number of 
bands (or particles). As in ID, the interesting part is 
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contributed by eigenstates at k^'s, because the eigen- 
state at k inv on each band is either parity even or parity 
odd. The contribution from these kj n „ 's to the Z 2 index 
is given by 

Dim(V 0Cc .) = n + (k mv ), (34) 

k in „ 

Dim(fy occ .) = ^jT n~(k inv ). 

k;„„ 

The Z 2 index obtained in a translationally invariant sys- 
tem remains unchanged in the disordered system as long 
as the disorder potential preserve inversion symmetry, as 
mentioned previously. We may generally state that the 
number of 1/2-in-gap states in the entanglement spec- 
trum for any inversion symmetric cut with or without 
weak disorder is given by the same Eq. (fT2"j) . with the 
only difference that now ki„„ runs over four and eight 
points in 2D and 3D respectively. 

In addition to giving the number of 1/2-in-gap states, 
the Z 2 index also gives the parity of Chern number in 
2D. In fact, we have 

= (_l)£>s„„ "~< k — ) (35) 
= (-1)°. 

In the last equality, we have applied Eq. (92,93) in 
Ref . [lH , which relate the Chern number to the inversion 
eigenvalues at all inversion invariant points in BZ. 

In 2D and 3D, if in addition to inversion symmetry, 
we also have time-reversal symmetry, every band is at 
least doubly degenerate (for spinful fermions) . And since 
[T, P] = 0, the two degenerate bands must have equal 
inversion eigenvalue at kj W 's. In this case (zi — 
must be an even number and one can further have 

(-1)(*1-*=V 4 = (-l)£»Nn. «~( k — )/ 2 (36) 

= (-ir°, 

where 70 is the Zi index of a TRI insulator. 

As a final remark, the Z 2 -index only concerns the to- 
tal number of occupied states with even/odd parities 
and therefore does not exhaust the topological classifi- 
cations. For example, an inversion symmetric insula- 
tor with only one filled band having B(T) = B(X) — 
—B(Y) = —B(M) = 1 is topologically different from an- 
other one with B(T) = B(M) = -B(X) = —B(Y) = 1, 
though they have the same Z 2 -indices. 

IV. ENTANGLEMENT SPECTRUM IN C n 
INVARIANT INSULATORS 

We are already aware that we can define a Z 2 index 
for inversion invariant insulators. This index represents 
the number of occupied eigenstates in each of the two 
ID representations of the inversion point group. In to- 
tal, there are n ID representations of point group C„, 
in each of which C n is represented by a root of the 



equation x n = (— 1) F . We denote the mth root by 

X m = e ^F+2(m-l))iv/n^ Tmg aUowg ug %Q & 

index [z\, z n ), in which z m is the number of occupied 
eigenstates in a subspace \&( x ">) 

z m = D*ro(*<£»>). (37) 

Any state \ip) in vf^™) satisfies C n \ip) = x m \?p}, or, in 
other words, transforms according to the m-th represen- 
tation of group C n . 

Yet it is not clear in what manner this Z n index is 
related to the number of protected in-gap states in the 
single particle spectrum. To answer this, we must de- 
fine what a symmetric cut is in a C n invariant insu- 
lator, analogous to the cut defined with inversion in- 
variant insulators, because a given entanglement spec- 
trum is associated with a specific real space cut. Due 
to the complexity of C n compared with the inversion 
point group, there are more than one type of symmet- 
ric cuts. Two integers mi and 7712, both being factors n, 
describe a symmetric cut denoted by A™, 2 mi that satisfies 
(a) the number of sites in A is exactly N/m\. (b) Sub- 
set A is invariant under m2-fold rotation, or symbolically, 
C m2 r £ A, Vr £ A and (c) the subset A, together with all 
copies of A obtained from acting C mim2 rotation on A for 
k < mi times, constitute the whole system, or symboli- 
cally, Ujfe = o,...,mi-iCTOima^ = ^' wri cre L represents the 
whole lattice. Notice that properties (a) and (c) together 
imply (d) C% ima A n C% im2 A = 0, if fa ^ k 2 mod mi, 
while C™ i 1 m2 = C m2 . In other words, an mim2-fold ro- 
tation on A generates other equivalent subsystems that 
are invariant under ?7i2-fold rotation. For n = 2,3,4,6, 
all possible combinations of (mi, m-i) are listed with typ- 
ical examples of every type of cuts in the real space in 
Table fl] Additionally, we note that for a fixed n, there 
is a one-to-one mapping between one type of symmetric 
cuts and one coset decomposition of a subgroup of C n . 

In the last column of Table [H the number of the pro- 
tected in-gap states is given in terms of Z n indices. 
This gives the least number of eigenvalues in the en- 
tanglement spectrum that nvv constrained in the range 
[1/mi, 1 — 1/mi]. If mi = 2, they are all at exactly 1/2 
as we have seen in the inversion invariant insulators. But 
when mi > 2, we do not have protected states at exactly 
1/2, but states that are prevented from moving outside 
the range of [1/mi, 1 — mi]. 

To prove these expressions, we first show that the en- 
tanglement matrix calculated with an A™ 2 mi , C{A™ 2 mi ), 
can be block-diagonalized into mi blocks. Since A is in- 
variant under C m2 , we can define a proper set of A with 
exactly Na/tui sites, A , such that C^r^Ao if r e A 
for p = 1, m 2 — 1. For every site r £ A we define the 
following rri2 symmetry adapted operators 




(y* q YC p m Mr)C-V 2l (38) 



where y q — cxp[i(2q + F — 2)7r/m 2 ] is the gth root of 
the m2 roots of e %Flx (q = 1, ...,7712). (For an intuitive 
schematic of this definition, see Figl3]) 
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FIG. 3. A picturesque expression of the definition of operators 
d 



Vq, 



We use A\j 2 as an example, and show how the opera- 
tors in A can be decomposed into three sectors with different 
eigenvalues of C3. 



These new operators have two 
{dy q , a (r),dl , ij9 (r')} = 5 qq >S a (s5 rr , and 



C'm 2 dy q ,a(Y)C' m2 — y q dy q _ a (r). 



properties: 



(39) 



The first property ensures that the original operators 
c a (r e A) can be linked to d y a (r <G Aq) by a unitary 
transform, thus preserving all eigenvalues of the correla- 
tion matrix. The second property leads to 



(40) 



To see this, we first notice that the ground state is C m2 - 
invariant, i.e., C m3 10) = e t6 \0), where 8 is some arbitrary 
angle. Then use Eq. (j3"9"j) : 



« 1 ,a( r X 2 .^ r ')) 
(0 1 C m l (Cm 2 dy n jCt (r) C m l ) (C, 



(41) 



n 2 dy q2 ,p(r')C m l)C m2 10) 



If qi 7^ 52, from the above equation, wc have 
(dt qi Jr)d yq2 Ar'))=0. 

From these we know that using the basis of d y |Q (r€ 
Ao), the correlation matrix is block diagonalized into 771.2 
blocks, i.e., 



C(A) = C (1 \A ) © C (2) (A )... © C (m2) (A ), (42) 



where 



C^ r ^{A ) = (dl ia (r)d Vp Ar'))- (43) 



Using this result, the study of in-gap states in C(A) re- 
duces to the study of in-gap states in each block (A ). 
When we have calculated the number of in-gap states, 
N^j d , for each block, we can simply add up all these 
numbers and obtain 



E 



N 



0) 



in id 



(44) 



N mid (AT/ mi ) ■ 

q=l,...,m 2 

for the total number of in-gap states. 

From here we only focus on one block, C^(A ), of 
the block diagonalized correlation matrix. In Appendix 
[Cl we prove that C^(Aq) can be written as a linear 
average of mi projectors 



C^(A ) = —(D^ 

777l 



£)(P)\ 



(45) 



where D r is a projector, projecting any state into a 

(v) 

subspace $f , in which every state satisfies these two 
conditions: (i) it is an eigenvector of C m2 of eigenvalue y p , 

and (ii) it is an eigenvector of C mim2 with eigenvalue . 
Notice that since C m2 = , we have the constraint 

(A[ p) ) mi = y p . Solve it and we have \ { r p) = exp[in{F + 

2((r— 1)tt72 +p— l))/(mi777,2)] for r = 1, m\. Given an 

(v) 

eigenstate of C n of eigenvalue x m =i,..., n , if it is in $r , 

using C mim2 = Cn 1 " 12 x m must satisfy Xm 1 ™ 2 = X'jfi. 

In Appendix [Cj we prove that the number of non-zero 
eigenvalues of is 



dim(Df>) = Dim{U Xm {W*™>\x, 

E 



(im)| T i>/(nimj) 



J2TT({r-l)m2+p-l+F/2)/(m 1 m 2 ) 



}) 



(46) 



^((r-llmj+t-l+F/Sl/lmjmj) 



Recall that in inversion invariant insulators, the differ- 
ence between dim(D) and dim(D) corresponds to the 
number of in-gap states; here we have a very similar re- 
lation. The block of the correlation matrix has 

N%l(A ) = . . max \dim(D^) - dim(D^)\ (47) 

cigenstatcs within the range [1/7711,1 — 1/777^ (see Ap- 
pendix [D] for proof). Therefore the total number of in- 
gap states is the sum of the number of in-gap states in 



each block: 

N mid (A?Lj = £ .max \dim(D^) - dim(Df)\. 

1 1 * — ' ij-l,...,mi J 
p=l,...,m 2 



Substituting Eq. (|4"6"|) . which relates dim(Dr') to Z n - 
index, into Eq. (|4"5)l . we obtain the expression of the total 
number of protected in-gap states for the symmetric cut 
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A") 2 as shown in the last column of Tabic (TJ 

The discussion above applies to systems with and with- 
out translational invariancc. In the case when transla- 
tional invariance is present, we desire to elucidate how 
the Z n index related to the eigenvalues of symmetry op- 
erators at high symmetry points. To understand the na- 
ture of this link, we calculate the Z n index given a single 
particle Hamiltonian in k-spacc, or in the homogeneous 
limit with periodic boundaries. To be exact, we want to 

(x ) 

count the dimension of the linear subspace Voce, where 
x q = cxp(i(_F + 2q — 2)n/n) for q = 1, First we 

write symmetry adapted operators in terms of annihila- 
tion operators the Bloch states, 7 m (k): 



-0x,,m(k) = 




( X * q yC? 7m (k)C-*>, (49) 



where B n > is the sewing matrix corresponding to the n'- 
fold rotation. This expression is non-zero if and only if 

(B n '(ko))mm = x q n ■ It can be physically understood 
as follows: since C n > = CrJ n , and (B n '(k )) mm is the 
eigenvalue of C n > for the state at ko, it must be the n/n'- 
th power of x q to be in the subspace of ^^ x i\ Further- 
more, in the BZ there are exactly n/n' points that are C n ' 
symmetric but not C n symmetric (n' < n) , and therefore 
there are n/n' orthogonal states, the linear combinations 
of which give the n/n' eigenstates of C n with eigenvalues 
satisfying x q l/n = (B„'(k )) mTO . 

It is straightforward to show C n 4'x q ,mO<)C~ 1 = 
x q ipe,m(k). Therefore if ip x m (k) ^ and m e occ, then 

state -0| g m (k)|O) belongs to subspace V c?\ contribut- 
ing +1 to z q . With Eq. ([5T)|) . the counting is therefore 
completely determined by the sewing matrices at these 
high symmetry points, or the number of each eigenvalue 
of C n ' at these points. The result for n — 2,3,4,6 are 
listed in Table [TXJ In Appendix [FJ we give an example of 

I 



where m is a real tunable parameter. This Hamiltonian 
is ^-invariant, C4(k)H(k x , k y )C^ 1 (k) = H(—k y ,k x ), 
with C<4(k) = exp(ia z TT / A)e lkl . Mark that C4-matrix 



where 7 m (k) is the annihilation operator of the Bloch 
state at k on the mth band. This definition is a lit- 
tle different from that for d Vqta in Eq. (|38|) . though they 
transform in similar ways under rotation: (i) d Vq ^ a are to 
block diagonalize the correlation matrix, while i\) x iin 's in 
general cannot do this; (ii) d y , Q (r)'s are not annihilation 
operators of single particle eigenstates, while tj) x TO (k)'s 
are generically annihilation operators of a single particle 
cigenstate; (hi) as we will see later, ip x m can be zero 
at points of symmetry while d y . Q is always nonzero. We 
notice that for a generic k that is not invariant under any 
subgroup of C n , {ip Xqi . m (k),^l^ m ,(k)} = 6q iq2 6 mm ,. 
Any one of these these k points on each band contributes 
+1 to every component of the Z n index. Again this con- 
tribution is ignored as it is completely determined by 
the total number of particles and total number of bands. 
But if ko is invariant under a subg roup Cfi* of Cni 
C n 'k = k , we have, applying Eq. (|B5I) to Eq.(j49]), 



(50) 



I 

deriving the row with n = 6, mi = 2, mi = 3 in Table \H\ 
following the general rules given above. 

Substituting the dimension counting into Eq. (|48[) . we 
have finally linked the number of in-gap states in a C n 
invariant insulator with a symmetric cut A™ 2 m ^ to the 

number of different eigenvalues of C m (m dividing n) at 
all high symmetry points. Finally, it should be noted 
that although we use the sewing matrix, defined only 
for a homogeneous system, to calculate the Z n index, 
the components of the index do not change under any 
adiabatic shift of Hamiltonian including addition of C n 
symmetric disorder potential. Therefore, as far as the 
disorder strength is small compared with the insulating 
gap, the number of in-gap states given in Table |H] is pro- 
tected from changing only by C n symmetry but not by 
translational symmetry, and is uniquely determined by 
the Z n index. 

In FigUl we give two specific examples of a C4 invariant 
2D insulator showing protected 1/2-in-gap states in the 
entanglement spectrum, both with symmetric cuts of the 
type Ay 2 . The system is described by a tight-binding 
model with two bands: 



(51) 

I 

is k-dependent since the rotation center is not chosen 
at a lattice site, but the center of a plaqucttc. Using 
C 2 (k x ,k y ) = Ci{—k y , k x )Ci{k x , k y ), we have C 2 (k) = 



^ m (U)= £ (x-;"/">„<(k )) mm )* ]T (x;rc^ m (k )c-p, 

v t=0,...,n'-l p=0,...,n/n'-l 



H(k x ,k y ) = (m - cos(k x ) - cos(k y ))a z + sm(k x )a x + sm(k y )a. 
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n 




m 2 


Typical 
examples 


Minimal number of midgap states in 
the entanglement spectrum 


2 


2 


1 


H B 


W ~z 2 \ 


3 


3 


1 


w w 


max |Z( — z;\ 

ij — 1.2,3 




2 


1 




|z x +z 3 -z 2 -z 4 | 


A 


2 


2 


M SB 


|Zl -Zsl + 1^2 -Z 4 | 




4 


1 




1,7 = 1,2,3,4 




2 


1 


\ i \ 1 


|z 1 +z 3 +z 5 -z 2 -z 4 -z 6 | 




2 


3 


®® 


|zi -Z4I + l z 2 -Zsl + |zs -Zel 


6 


3 


1 




maxflz! + z 4 - z 2 - z 5 1, |z x + z 4 - z 3 - z 6 1, |z 2 




3 


2 




maxdz! - z 3 1, |Zi - z 5 1, |z 3 - z 5 1) 

+ max(|z 2 - z 4 |, |z 2 - z 6 |, |z 4 - z 6 |) 




6 


1 




max Iz,- — z,| 

i,/=l,...6 y 



TABLE I. All possible symmetric cuts A^? mi in C n invariant insulators. In the fourth column, examples of each case are shown 
and in the last column, the minimal numbers of in-gap states in the entanglement spectrum are represented in terms of the 
Z n -index. 



it7 z e ltykx ~ ky \ We assume that the system is half- filled. 
There are two ^-symmetric points T and M, where 
the Hamiltonians are H(T) = (m — 2)a z and H{M) = 
(m + 2)tr z , respectively. If m > 2, the occupied state 
is (0,1) T at both points, giving B{T) = e" l7r/4 and 
B(M) = -e~™ /4 . At X or Y, H(X/Y) = ma z , and the 
occupied state is still (0, 1) T , giving B(X/Y) — i. Ac- 
cording to Table ITU we find z\ = z% = Z3 = za = 1. For 
?77, < —2, similarly, we can similarly derive .21,2,3,4 = 1. 
For < m < 2, at r, the occupied state becomes 
(1,0) T , contributing B(T) = e l7r / 4 , while for M and 
X, the occupied state is still (0, 1) T . In this case we 
have z\ = 2, zi = 2:3 = 1 and z^ = 0. Finally for 
—2 < m < 0, one can similarly derive z\ — z± = 1, 
Z2 = 2 and 23 = 0. Summarizing all possibilities, we 
have Z\ = 1 — sgn(m — 2)/2 + sgn{m), 22 = 1 + sgn(m + 
2)/2 — sgn(m)/2, 23 = 1 — sgn{m + 2) + sgn(m)/2 and 
24 = 1 + sgn(m — 2)/2 — sgn(m). From Table UJ we 
know that with this type of symmetric cut, the number 
of 1/2 in-gap states are N 1 / 2 = \zi — 2 3 | + \z 2 — zs\ = 
2\sgn(m + 2) — sgn(m — 2)|. If |m| > 2, then there 
is no protected in-gap state and if |tt7| < 2, there are 
always two degenerate entanglement eigenvalues at ex- 
actly 1/2. In FigUl we confirm this result by calculating 
the entanglement spectrum with or without C4 invariant 
impurities. In the top two figures of Fig|3J we can see 
that although the details of the spectrum depends on the 
specific choice of the real space cut, the presence of the 



doubly degenerate protected in-gap states is common to 
both cuts. In the bottom two figures of FigHJ we can see 
that random impurities in general push the eigenvalues to 
the two ends of the spectrum as the electronic states be- 
come more localized, the doubly degenerate in-gap states 
always remain at 1/2. 

V. DISCUSSION 

A. Entanglement spectrum of 3D insulators with 
PGS 

The major focus of the paper is 2D insulators with 
77-fold rotation symmetry; however, the results obtained 
can be promptly extended to 3D insulators with 77-fold 
rotation symmetry. This is because a 3D C n insulator 
can be considered as a 2D C n insulator parameterized by 
a continuous k z along the rotation axis. For every k z we 
can find its number of in-gap states in the entanglement 
spectrum and, as we require the 3D system to be insulat- 
ing, this number of in-gap states must be a constant for 
different k z . Therefore, studying the system at one k z - 
slice as a 2D insulator is sufficient to understanding the 
symmetry protected in-gap states in the entanglement 
spectrum in the 3D C n -invariant system. 

For 3D insulators with other than C n PGS, their entan- 
glement spectra in principle show other symmetry pro- 
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n = 2 


Dim(W occ 


4 00 + 4 J (A)+4 >(Y) + n { 2 >(M) 




_. , T ( e t(*'+JWJ) 


(e i(*'+W3), {e i(F + 2),/2 ) , i(F + 2),r/2 ) {( ,i(F + 2)^/2 ) 

4 '00 + 4 J P0 + 4 '00 + 4 ( M ) 


n — 3 


Dim(Wcc ') 


4 > (r) + ' (A ) + J (A") 






n 3 e > (T) + n ( 3 ' (A) + n ( 3 ' (A ) 




£>im(tf£c j ) 


4 '00+4 , s O + 4 


n = 4 




r^ e '(r) + n^ e j (M) + n^ e > (X) 






nf > (r) + >{M)+ n\ e > (X) 






< (r)+< ; ( A 0+4 






00+4 j (M)+ni, 


n — 6 


Dimple ') 


p ' ] ✓ i \ IF ' 1 / t r-\ IP ' l/n*-\ 

4 00 + 4 i K ) + 4 ( M ) 






n£ '(r) + n^ ; (A) + 4 J (M) 






4 00 + 4 '0*0 + 4 W 






4 '00 + 4 '(lO+na '(AO 






/„i(F + 8)7r/6^ , ;(F + 2)tt/3n ,„iF7r/2s 

4 '(n+n^ 1 ' '>(*) + 4 (Af) 




«m(*i E c r +iu " /H) ) 


/ i(i'+10K/6l ( e i(F + 4)^/3> ( i(F + 2)^/2> 



TABLE II. Z"-index in a C„ invariant insulator with translational symmetry. These indices are related to the number of all 
eigenvalues of the corresponding sewing matrices associated with subgroups of C n at high symmetry points within the BZ. n^, 
represents the number of states with eigenvalue e of m-fold rotation C m , where C m is a subgroup of C n . 



tected features, and the number of in-gap states depend 
on other factors. Nevertheless, if the underlying PGS 
has only ID representations, as C„ group does, one can 
still define the Z 9 ° -index, where go is the order of the 
group, similar to the Z"-index for C„ -invariant insula- 
tors. Then one can make exact statements about the 
number of in-gap states for properly defined symmetric 
real space cuts. If the group has higher-dimensional ir- 
reducible representations, it is yet to be studied how the 
occupation numbers of states belonging to these repre- 
sentation are related to the entanglement spectrum. 



B. Effects of weak interaction 

Throughout the paper we have been studying systems 
without interactions. Now we discuss how weak interac- 
tion changes our results. By 'weak', we mean that the 
ground state of the interacting system can be obtained by 
adiabatically transforming the Slater determinant state 
or the ground state of the non-interacting system. This 
implies that the interacting ground state is both gapped 
and non-degenerate. Therefore, we are not considering 
any symmetry breaking state as such a state must be de- 
generate with at least another state related by the bro- 
ken symmetry. Additionally, we are not considering frac- 
tional quantum Hall systems the ground states of which 
are degenerate. In this work, we assume the interacting 
ground state to be an insulator with a non-degenerate 
ground state and a many-body gap. 

In an interacting insulator, single particle states lose 
their meaning, but the many-body ground state is still a 
ID representation of C ra , giving a Z n index. The relation 
between this Z n index and the Z n indices can be found by 



noticing that the C n rotation operator can be expressed 
as 

C n = exp[i^ ( m ~ 1 )4 m , Q ( r )^ m ,a(r)] 

m-l,...,n,Q,r£L/n 

(52) 

where 

d Xm ,a(r) = -7= ^ P m C p n c a (r)C-P, (53) 

p—l,...,n 

and that 

z m = ( ]T 4 m »<4 ra ,a(r)>. (54) 

a,r£L/n 

Since the ground state is an eigenstate of C„, we have 

2tt 

{C n ) = cxp(i — z) (55) 
n 

= cxp(i— VVm - l)z m ), 
n t— 1 

m 

or 

z = (m — l)z rn mod n. (56) 

m— l,...,n 

Eq. (|56"|) shows how the Z n classification in a non- 
interacting system downgrades to the Z n classification in 
an interacting one. This indicates that a non-interacting 
insulator with a non-trivial Z n index may have a trivial 
Z n index in presence of weak interaction, and hence two 
topologically distinct non-interacting insulators may be- 
come indistinct as interaction is turned on, similar to the 
situation in superconductors^ - — . 



14 



Homogeneous 



Disordered 



1 

0.8 

0.6 

0.4 

0.2 


1 

0.8 
0.6 
0.4 
0.2 




1 

0.8 

0.6 

0.4 

0.2 


1 

0.8 
0.6 
0.4 
0.2 




FIG. 4. Entanglement spectra in a Ca invariant insulator described by the Hamiltonian in Eq, (|51|l with m = 1.5 are plotted 
with two different cuts of the type A\/ 2 on an 8-by-8 square lattice. The vertical axis is the entanglement eigenvalue and each 
level in the spectrum is represented by a horizontal line and Different columns are for different specific cuts; the first row shows 
the result without impurity (with translational symmetry) and the second with random disorders of strength^ 0.4 that are C4 
symmetric. Red lines mark the position of eigenvalues at a value of exact 1/2. The number of protected in-gap states is two 
in all cases. 



Beyond the basic discussion presented here, the effects 
of interactions in point group symmetric topological in- 
sulators have many unanswered questions. A particu- 
larly important question is how the effects of interactions 
present in the changes in the entanglement spectrum. It 
is not clear how the in-gap states in the entanglement 
spectrum evolve as one turns on the interaction. We 
leave this question and other questions as to the role of 
interactions to future work. 



VI. CONCLUSION 

In our analysis of the resultant entanglement spectrum 
in point group symmetric topological insulators, we have 
proven the relation between the 1/2-in-gap states in the 



entanglement spectrum in an inversion invariant insula- 
tor and the inversion eigenvalues at inversion invariant 
points. This was accomplished by noticing that, given 
any inversion invariant non-interacting insulator, one can 
always find a complete set of single particle eigenstates, 
of both the single particle Hamiltonian and inversion op- 
erator, and the numbers of occupied eigenstates that are 
odd/even under inversion are good quantum numbers. 
This results in a Z 2 -indcx of the insulator. The sum 
of these two integers gives the total number of particles 
and their difference is exactly the number of 1/2-in-gap 
states in the entanglement spectrum. We have extended 
this result from inversion invariant insulators to 2D in- 
variant insulators with C n invariance: The number of 
occupied states that transform according to each ID rep- 
resentation of C„ is a good quantum number, resulting 
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in a Z"-index. In a translationally invariant system, this 
Z n -mdex can be calculated by counting the number of 
each C m (m dividing n) eigenvalue at high symmetry 
points. In the entanglement spectrum, the number of 
protected in-gap states is shown completely determined 
by the Z n index with or without translational invariancc 
in any C n invariant insulator. Specially, if and only if all 
components in the Z n index are equal, then the resul- 
tant number of protected in-gap states is zero. Finally, 
we briefly discussed how a similar analysis can be ex- 
tended to 3D insulators with point groups having only 
ID representations, and demonstrated that under weak 
interaction, the Z"-index for a C n invariant insulator re- 



duces to a Z n -indcx. 
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Appendix A: Subtleties in the definition of inversion operator in k-space 

The relation between the number of in-gap states and the parities of occupied states at k = 0,7r in a ID inversion 
invariant insulator, Eq. (|T2| . is completely consistent with Eq.(12) of RefHH, although two expressions take different 
forms. 

The difference lies in the definition of the inversion operator. While RefJH uses a fc-independent V as the matrix 
representation of inversion operator in fc-space, in this paper we use a fc-dependent V(k) = Ve lk as in Ref.JHj]. The 
factor e lk comes from the fact that the inversion center is mid-bond, instead of on-site. Therefore, at k = n, the 
even/odd parity states in our paper correspond to odd/even states in Refl32l. 

In 2D, the inversion operator in k-space is P(k) = r Pe lk:c+ ' lkv on a rectangular lattice, because the inversion center 
is at (ai +a2)/2; in 3D, we use V(k) = Ve tfe *+ lfe v+ tfe * as the inversion matrix in k-space, because the inversion center 
is set at (ai + a.2 + s^)/2. An advantage of using such definitions is that the formula for the number of in-gap states 
in all dimensions takes exactly the same form: 

Nl/2 = | Tl +( k t"«) _ n -( k mv)\- (AI) 

kirn, 

Appendix B: A basic property of the sewing matrix 

In this Appendix we derive a simple property of a general sewing matrix associated with a PGS operation R we 
have applied to prove Ea.([25|). The sewing matrix Bft(k) is defined as 

B mn (k) = (u m (Rk)\R\u n (k)}, (Bl) 

where m,n £ occ. Multiply both sides by |u m (i?k)) then sum over m G occ, we have 

B mn \u m (Rk)) = \u m (Rk))(u m (Rk)\R\u n {k)). (B2) 

Since R is a symmetry of the Hamiltonian, if one i?|u n (k)) must also be an occupied state at Rk if n G occ, therefore 
J2m£occ \ u m(Rk)) (u m (Rk)\ can be replaced by the identity in the occupied bands. Eq. (jB2[) becomes 

fl|u n (k))=B mn |u ro (.Rk)). (B3) 

Then we substitute 

|n„(k))= 7 t(k)|0), (B4) 
\u m (Rk)) = jl(Rk)\0) 

into Eq. (|B3p . and obtain 

R^nWR- 1 = B mnlm {Rk). (B5) 
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Appendix C: Decomposition of a correlation matrix into projectors 

In Sec. IIV1 we have shown that the correlation matrix of a A™?^ cut can be block diagonalized into m.2 blocks, 

each denoted by C^(A). In this appendix, we prove that this matrix can be further decomposed into the average of 
mi projectors: 

C^(A) = — (£>< s) + ...+£>W). (CI) 
mi 

(s) 

We further prove that the number of unity eigenvalues in the projector D\ equals the number of occupied states the 
C n eigenvalues A's of which satisfy \ n /( m ^) = e a((i-i)m J +»-i+F/2)/(mim J ) i 

Let us start from recalling that c\ s J ^ = {d^ yaja {ri)d Vs ^{rj)), where r^j G Aq. For each d Vs ^(vi), we further define 
mi new operators 

/i, s 2(r I ) = -L E ( e *) P ^ 1 ™ 2 4 s V 1 )^ im2 P , (C2) 

V p-0,...,mi-l 

where e mi = exp[z 2tt( s - 1+F/2 )]. There are in total mi e's that satisfy this condition, which are denoted by e& = 
e i2ir(km2+s+F/2)/(m 1 m 2 ) f or _ ..^ mi, and are all eigenvalues of operator C mi m 2 - This is very similar to what we 
have done in the ID inversion symmetric insulators (see Sec[njA)): over there d a (r) and d a (r) are combinations of 
c Q (r) that are even and odd under inversion; here /e,2(rj) are linear combinations of d y ^. a (r) with different eigenvalues 
of C mim2 . These operators have two properties (a) C , TOim3 /i,a(ri)C'+, ima = e/ififc) and (b) {/i.afo), / w e , ^(r,-)} = 
See'&apSnTj- From the first property, one derives (/^ e Q ( r i)/i' S ^( r j)) = if e 7^ e'. The derivation is similar to that 
of Eq. pU)) . only that here we use the fact that the ground state is an eigenstate of C mim2 . Therefore the correlation 
matrix C^ s \Aq) decomposes in the form of Eq. (|45|) where 

(i^W/J = </Wl ka (rO/W (r 3 -)). (C3) 

For simplicity, the superfix (s) will be suppressed and always implied throughout the remainder of this Appendix. 

Now we need to prove that for every k = 1, ...,mi, matrix Dk is a projector, i.e., (Dk) 2 = Dk- First we notice 
that Ctti 1 7712 is a subgroup of C n or C„ itself, so each single particle eigenstate has an eigenvalue of C mim2 . There are 
m\m,2 eigenvalues of C mim . 21 dividing the single particle Hilbert space into m\m,2 sectors, each of which is denoted 
by \l/ e = \&£ cc + ^u nocc , where e is an eigenvalue of C mim2 . Suppose that in the e-sector, all eigenstates are generated 
by tpl q , with q = 1,...,N * N orb /(m 1 m 2 ). The operators V'l.g's tnus satisfy (a)C mim2 ipl q C} nim2 = etpl q and (b) 

The previously defined operators f €k! a( r i) can be decomposed into a linear superposition of all single particle 
eigenstates: 

/ £fc ,a(rO - E W$ q i> x , q , r, G Ao. (C4) 

A 

But according to the property (a) of /e fc ,a(ri) and property (a) of ij)\,q, we have ^ only if A = and therefore 

A k ,a(ri) = W&e,,*. (C5) 

One is reminded that rv G Ao, which is l/m.2 of A and is therefore l/(mim2) of the whole lattice, therefore W^ e ' 
is in general a N * A or f,/(7nim2)-dimcnsional unitary matrix. Before proceeding, one is reminded of the difference 
between the matrix U introduced in SecUHA), which is a unitary transform between the unsymmetrized c a (r)'s and 
the symmetrized single particle eigenstate operators, and the matrix W, which is a transform between the symmetrized 
operators f ejOI and the symmetrized single particle eigenstate operators. In general, W is a block-diagonal matrix with 
each block labeled by its A. Each block in a block-diagonalizcd matrix must be unitary, if the whole matrix is unitary. 
Using Eq. (|C5[) . one can calculate Dk'- 

{D k ) ia , Jl3 =£^K^K7W<C ? VW> (C6) 

q,q' 

= V iy (efe)t ■ W^hn a 

/ j q,ia" JP,q- 

q£Locc 
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And D\: 

{Dl) ia jj3 = ^(Ufc)ia, m7 (£>fc)m 7 j/3 (C7) 

— 2_j vv q,ia vv mf,qvy q',my yy JP,q- 
q 7 q' £occ 

Using the unitarity of we have 

W^ m ^ q W^\,, mi = (W^W^) qq , = S qq ,. (C8) 
Substituting Eq.JC8} into Eq.|C2|, we have 

(D 2 k ) ia , j(i = Y, W^l ia W^ jP , g = (D k ) ia>j0 . (C9) 

q£occ 

From the explicit expression of Eq. (|C6[) . it is straightforward to see that for every q € occ, the column vector 
^ ia q ' 1S an eigenvector of of eigenvalue 1 and for every q ^ occ, it is an eigenvector of eigenvalue zero. 
Therefore dim{D k ) is equal to the number of occupied states with C mim . 2 eigenvalue e^, or the number of occupied 
states with C„ eigenvalue A that satisfies A m i m = = e^. 



Appendix D: Number of in-gap eigenvalues of a matrix C = —(D\ + ... + D n ), where Di is a projector 

In this Appendix, we prove that if a matrix C can be written as the average of n projectors, i.e., C = \jn{D\ + 
... + D n ) 1 then there must be at least N m id eigenstates, the eigenvalues of which, e, are in the range 

ee[l/n,l-l/n], (Dl) 

. The number of these states is given by 

N mid = max \Dim{Di) - Dim(Dj)\, (D2) 

i,j=l,...,n 

where Dim(Di) is the number of nonzero eigenvalues of Di. 

Before going on, the authors would like to point out that the proof is completely mathematical and has very little 
relation to the physics presented in the paper. It is provided here only for completeness of the work. 

First we prove a lemma: Given two hermitian matrices A and B, with eigenvalues a^s and b^s, consider their 
sum M = A + B, then for any eigenvalue of A, say a^, there must be an eigenvalue of M that satisfies + b m i n < 
m* < at + b max , where b m i n (b max ) is the minimum (maximum) eigenvalue of matrix B. We can always define 
B' = B — (b m i n + b max )/2I 1 while B and B' have the same eigenstates, and here the maximum/minimum eigenvalue 
of B' is ±(b ma x — bmin) /2 = ±&o- So for any eigenstate of M' = A + B' whose eigenvalue is between a; — bo and 
cii + bo, the same state must also be an eigenstate of M = A + B with eigenvalue between a; + b m i n and + b max . 

Denote the eigenstate of A with eigenvalue a with |a), then consider the quantity 

(a\B' 2 \a)=J2\y 3 \ 2 b' J 2 <bl (D3) 

where we have used the expansion of B' in terms of its eigenstates B l = ■ bj\bj)(bj\ and the expansion of \a) in \bj), 

On the other hand, one can always decompose |a^) in terms of eigenvectors of M 1 

\a,) =Y x n\m' n ), (D4) 

71 

where \m' n ) is an eigenvector of M' with eigenvalue m' n . We again calculate the quantity 

(a\B' 2 \a) = (a\(M' - A) 2 \a) (D5) 
= (a(M' -af\a) 
= YW\ 2 {m' n -af 

n 

<b 2 . 
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Since \x n \ 2 > for each n and J2 n \ x n\ 2 = 1, there must be at least one n* that satisfies 

\m! n , -a\< b , (D6) 

or 

a — bo < m n * < a + bo. (D7) 

Therefore for we know there must be an eigenvalue of M, m n * = m' n , + (b max + b m i n )/2 that is between aj + b m i n 
and at + bmax. 

Consider the matrix D = D\ + ... + D n , define A = D\ + D 2 and B = D3 + ... + D n . If dim(Di) ^ dim(D2), 
then there must be at least £12 = \dim(D\) — dim(D2)\ number of eigenstates of A with eigenvalue 1. On the other 
hand, since each Di is a projector, b m i n = and b max = n — 2. Suppose \(f>i,...,s 12 ) are the degenerate eigenstates 
of A with eigenvalue 1, and they make a subspace called $>a- Take |ai) = |</>i), from the lemma, we can prove that 
there must be one eigenstate |Ai) with eigenvalue between [l,n — 1]. Then if 612 > 1, one is guaranteed to have a 
state within ^a that is orthogonal to |Ai). Take that state as \a,2}. The decomposition of \a2) does not contain |Ai), 
or, x\ t = 0. Therefore, for Eq. (|D5|) to hold, there must be another eigenstate of M, called | A2) that has eigenvalue 
between [1, n— 1]. Then if 612 > 2, one is guaranteed to have a state within $^ that is orthogonal to both lAi^), called 
1 03). Start from |a3), one obtains another eigenstate of M with eigenvalue between [l,n — 1]. The process continues 
until one has exactly S12 eigenstates of M with eigenvalues [l,n — 1]. Now, one is not guaranteed to have a state in 
$A that is orthogonal to all | Ai^. . . ,<s 12 ) - In this way, we have proved that there must be at least \dim(D{) — dim{D2)\ 
eigenstates of D having eigenvalues between 1 and n — 1. It means there are \dim{D{) — dim{D2)\ number of protected 
in-gap state in the spectrum of C = D / 'n. 

Notice that in the proof, D\ and D2 are arbitrarily chosen, and we can in principle choose any two D^s. This 
means the least number of protected in-gap states is 

Nmid = max \dim{Di) — dim{Dj)\ 1 (D8) 

i,j—l,...,n 

and these states are protected because their eigenvalues must stay within the range [1/n, 1 — 1/n] and can never 
approach or 1 infinitely. 



Appendix E: Explicit derivation of the number of in-gap states for a A^ 2 -cut in a Ce-invariant system 

In this Appendix we explicitly show how we derive the expression of the number of protected in-gap states in terms 
of the Z n index for a A\^ 2 cut in a Cg-invariant insulator. 

In this system, the index z m =i 1 ... 1 6 represents the number of occupied states in the m-th representation of C§, that 
is, the states that have eigenvalue exp(i7r(F + 2m — 2)/n) under C§. According to the process sketched in Sec lIVl first 
we can block diagonalized the correlation matrix C(A^ 2 ) into three blocks. Then the total number of in-gap states is 
simply the sum of the numbers of in-gap states in each block. In the first block describes the entanglement between 
states with C3 eigenvalue exp(iFir/3) inside and outside A; the second block describes the entanglement between 
states with C3 eigenvalue exp(in(F + 2)/3) inside and outside A; the third block describes the entanglement between 
states with C3 eigenvalue exp(i7r(F + 4)/3) inside and outside A. 

Next, we notice that states with the same C3-eigenvalue, i.e., contributing to the same block, in general have two 
different C 6 -cigenvalues. States have C 3 -eigenvalue ex-p(iFir/3) may have C 6 -eigenvalue exp(iF7r/6) or exp(i(F + 
6)71-/6); states with C 3 -eigenvalue exp(i(F + 2)n/3) may have C 6 -eigcnvalue exp(i(F + 2)tt/6) or exp(i(F + 8)71-/6); 
states with C3-eigenvalue exp(i(F + 4)7r/3) may have Cg -eigenvalue exp(i(F + 4)7r/6) or exp(i(F + 10)7r/6). Therefore, 
an occupied subspace with a certain C3-eigcnvalue can be split into two subspaces with different Cg-eigcnvalues. The 
number of in-gap states in each block is exactly equal to the difference between the number of states in these two 
C6-subspaces. 

According to this analysis, we have N^j d = \z\ — z 4 |, N^J d = \z 2 — z 5 \ and N^J d = \z$ — Zq\. The total number of 
in-gap states is N mid = N^ d + N%} d + N^ d = \z x - z A \ + \z 2 - z 5 \ + \z 3 - z 6 \. 



Appendix F: Explicit derivation of Dim^lcc + ' ) in Table HT1 

if f + 4)tv /G 

In this Appendix, we use the example of counting the dimension of ^% cc to illustrate the way we use to derive 

Table HI 
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From Appendix \C\ we know that Dim{^% cc + 8 ) is the number of occupied states whose C6-eigenvalue is 
e t(F+4)7r/6 -^/g hayg a i so shown that only high symmetry points contribute to the difference of the dimensions 
of subspaces, and therefore we ignore the contribution from states at generic k's. 

At r, the highest symmetry is Cg, so the contribution to Dim( i i!^ +4)7 ' /a ) is simply the number of states with 
Cq eigenvalue e I ( F + 4 ) 7r / 6 . At K, the highest symmetry is C3. Notice that C3 = C|, and we can see that if the C3- 
eigenvalue of a state at K is e J ( F + 4 W 3 j then it will contribute +1 to Dim^l^^ 4 ^ /e ) . Finally, M is a ^-invariant 

point. Notice that C2 = Cf , and we can conclude that a state at M contributes +1 to Dim(^^ c + 1 ) if and only 
if its ^-eigenvalue is cxp(iF7r/2). 
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